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In addition to gamma rays, dark matter annihilation products can include energetic electrons 
which inverse Compton scatter with the cosmic microwave background to produce a diffuse ex- 
tragalactic background of gamma rays and x rays. In models in which the dark matter particles 
annihilate primarily to electrons or muons, the measurements of EGRET and COMPTEL can pro- 
vide significant constraints on the annihilation cross section. The Fermi gamma ray space telescope 
will likely provide an even more stringent test of such scenarios. 
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INTRODUCTION 

Annihilations of dark matter particles are predicted to 
produce a wide variety of potentially observable prod- 
ucts [l[. These include neutrinos, antineutrinos, elec- 
trons, positrons, protons, antiprotons, and photons. 
Searches for the photons generated in dark matter anni- 
hilations have a major advantage over most other indirect 
detection techniques in that these particles travel essen- 
tially unimpeded. In particular, unlike charged particles, 
photons are not deflected by magnetic fields, and thus can 
potentially provide valuable angular information. Fur- 
thermore, below approximately ~10 GeV, gamma rays 
are not attenuated, and can propagate without energy 
losses over cosmological distances. 

The transparency of the universe to gamma rays opens 
the possibility of observing a difl^use, isotropic back- 
ground associated with dark matter annihilations from 
throughout the observable universe. Dark matter parti- 
cles with weak-scale masses which annihilate to quarks 
or gauge bosons lead to a diffuse spectrum of energetic 
photons which peaks at ~1-10 GeV and could poten- 
tially be studied with the Fermi gamma ray space tele- 
scope (FGST) [aK This possibility has been studied else- 
where in detail [3, 01 . Recent observations of the cosmic 
ray positron fraction by PAMELA [5], and the cosmic 
ray electron (plus positron) spectrum by ATIG @ and 
FGST , have generated a great deal of interest in dark 
matter particles which annihilate primarily to leptons 
[loj . The energetic electrons produced in such annihila- 
tions rapidly transfer their energy via inverse Compton 
scattering into lower energy photons, such as those con- 
stituting the cosmic microwave background. The final 
product of this process is a diffuse extragalactic back- 
ground of lower energy gamma rays and x rays. 

In this paper, we study the extragalactic diffuse 
gamma ray background resulting from the inverse Comp- 
ton scattering of energetic electrons produced in dark 
matter annihilations. We find that if the annihilation 



cross section (or equivalently, the boost factor) is chosen 
to generate the signal observed by PAMELA, ATIC, or 
FGST, a very bright isotropic gamma ray background is 
predicted. Existing measurements from COMPTEL 11 1 
and EGRET [13] can be used to constrain the dark mat- 
ter annihilation cross section. Observations by FGST 
will provide a powerful test of dark matter annihilating 
to electrons or muons. 



THE EXTRAGALACTIC DIFFUSE GAMMA 
RAY BACKGROUND FROM DARK MATTER 
ANNIHILATIONS 

The dark matter annihilation rate per volume at a red- 
shift z is given by 



Riz) 



Mn 



— (M,z)il+zrdMy 



X 



p^{r,M)A'Kr'^dr, 
(1) 



where dn/dM is the differential comoving number den- 
sity of dark matter halos of mass M, {av) is the dark 
matter annihilation cross section, mx is the mass of the 
dark matter particle, and p is the density of dark matter 
in a halo as a function of the distance from the center of 
the halo, r. 

The number density of dark matter halos of a given 
mass as a function of redshift is given by 
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where pM is the average matter density, <t{M, z) is the 
variance of the linear density field, and f{<J~^) is the 
multiplicity function. The redshift and cosmology de- 
pendence is contained in cr{M, z), which can be defined in 
terms of the matter power spectrum, P{k), and top-hap 
function, W{k,M) = {'i/k^R^)[siTL{kR) - kRcos{kR)], 
where R = (3A^/47rp„J)^/^ 
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In determining the cold dark matter power spectrum [13| , 
we adopt flbh'^ = 0.02267, Qch'^ = 0.1131, = 0.726, 



and h = 0.705, as measured by WMAP [lj|. The growth 
function, D{z), is the linear theory growth factor, nor- 
malized to unity at z = [l^. We use the ellipsoidal 
(Sheth-Tormen) form of the multiplicity function j 
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where p = 0.3, Ssc = 1-686 and a = 0.75 11]. We fix A = 
0.3222 by the requirement that all of the mass resides 
in halos. The halo mass function is most sensitive to 
variations in the cosmological parameters erg and Ug, for 
which we adopt ris = 0.96 and as = 0.812. 

To describe the distribution of dark matter particles 
within halos, we use the Navarro-Frenk- White (NEW) 
profile 18], and halo concentrations given by the analytic 
model of Bullock et al. [l^ . Rather than extrapolating to 
lower masses, however, we fix the halo concentrations to 
a constant value below M = 10^ Mq (as in Ref. Dif- 
ferent choices for the halo profile or concentration could 
plausibly increase or decrease the annihilation rate and 
gamma-ray flux, but likely only by a factor of a few or 
less. 

Following Ref. ]^4j, we can write the flux of photons 
from dark matter annihilations throughout the observ- 
able universe in the form 



scattering. Since background radiation evolves with red- 
shift, the spectrum of inverse Compton photons is cal- 
culated as a function of redshift and added to the spec- 
trum of prompt photons. To calculate the prompt pho- 
ton spectrum from dark matter annihilations, as well as 
the prompt electron/positron spectrum, we have used 
the program PYTHIA [2l[. For details of the inverse 
Compton spectrum calculation, we refer the reader to 
the Appendix. The results in the next section include 
the spectrum of inverse Compton photons along with the 
spectrum of prompt photons. 



RESULTS 

In Fig. [U we show the spectrum of the extragalactic 
diffuse gamma-ray background from dark matter anni- 
hilating to each species of charged lepton (e+e^, /i^/i^, 
T+r^). For comparison, we also show the results for a 
dark matter particle annihilating to W~^W~. In each 
frame, the curves denote dark matter masses of 100, 200, 
400, 800, and 1600 GeV, from left to right. In each 
case, we have assumed an annihilation cross section of 
(av) = 3 X 10~^^ cm^/s, which is the typical value for 
a thermal relic. Also shown in each frame are the mea- 
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where px denotes the average density of dark matter, 
A^(z) is the average squared overdensity and h{z) = 
+ rim(l + z)^ with f^A and 11m being respectively 
the present fractions of critical density given by dark en- 
ergy and matter. The function r describes the estimated 
optical depth of the universe to gamma-rays [2^ . 

The spectrum of gamma-rays per annihilation, 
dNj/dEj, depends on the dominant annihilation chan- 
nels, and on the annihilating particles' mass. Neutrali- 
nos, for example, typically annihilate to final states con- 
sisting of heavy quarks or gauge/Higgs bosons. While the 
decays of such particles produce a spectrum of prompt 
photons, this spectrum is rather soft. In contrast, the 
gamma-ray spectrum from a dark matter particle anni- 
hilating to leptons can be significantly harder {i.e. de- 
positing a larger fraction of the energy into photons with 
energies not far below the dark matter particles mass). 
For the additional contribution from final state radiation, 
see Ref. (2^ . 

In addition to prompt gamma-rays, dark matter an- 
nihilations can also produce high energy electrons and 
positrons. Such electrons/positrons transfer their energy 
to the background radiation through inverse Compton 



surements from the COMPTEL P| and EGRET [T2| 
experiments. 

In Fig. [21 we again plot the diffuse gamma-ray back- 
ground from dark matter, but have increased the an- 
nihilation cross section (or boost factor) to the maxi- 
mum value consistent with the measurements of COMP- 
TEL and EGRET. For dark matter particles that anni- 
hilate to W~^W~ (or other varieties of heavy fermions or 
gauge/Higgs bosons), the strongest constraints typically 
come from the prompt gamma-ray emission, in contrast 
to the inverse Compton photons. A 100 GeV (1 TeV) 
dark matter particle annihilating to W^W~ would ex- 
ceed egret's diffuse background measurement if the 
annihilation cross section were boosted more than ^^400 
7000) beyond the common thermal estimate (3 x 10~^^ 
cm^/s). Dark matter that annihilates to leptons (espe- 
cially e+e~ or /i+/x~) deposit a larger fraction of their 
energy into electrons, and thus ultimately into lower en- 
ergy inverse Compton photons. In the upper left frame of 
Fig. [21 we find that we cannot increase the annihilation 
cross section of a 100-500 GeV dark matter particle an- 
nihilating to e+e^ by more than ~200-500 relative to the 
thermal estimate without exceeding the measurements of 
COMPTEL or EGRET. In Fig. [1 we show these maxi- 
mum allowed boost factors for each annihilation channel 
as a function of the dark matter particle's mass. The 
inverse Compton photons described here may also con- 
tribute significantly to the reionization of hydrogen and 
helium gas at redshifts of z ^ 6 — 20 24] . 

These results are very interesting within the context 
of recent observations from PAMELA 5] , ATIC @ , and 
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FIG. 1: The diffuse extragalactic gamma-ray background from dark matter particles annihilating to e e~, ^ r+r~, or 
W^W~ . In each frame, the curves denote dark matter masses of 100, 200, 400, 800, and 1600 GeV, from left to right. In each 
case, we have assumed an annihilation cross section of {av) = 3 x 10~^® cw? /s (and no boost factor), which is the typical value 



for a thermal relic. Also shown in each frame are the measurements from the COMPTEL [Till and EGRET [Is 
(on the left and right, respectively). 



experiments 



FGST 0] , each of which have been interpreted as possible 
indications of dark matter annihilations taking place in 
the Galactic halo [1, 13, 25 1. In order for these signals to 
be the product of dark matter annihilations, however, the 
annihilations must proceed largely to charged leptons, 
and at a considerably higher rate than would be naively 
predicted for a typical thermal relic distributed smoothly 
throughout the Milky Way. Thus if one interprets these 
signals as a product of dark matter annihilation, then a 
significant diffuse background of gamma rays from the 
inverse Compton scattering of dark matter annihilation 
products is also predicted. 



As a first example, consider a ^ 300 GeV dark matter 
particle which generates the PAMELA positron excess 
by annihilating directly to e^e^ with a boost factor of 
~30-250 In this case, the inverse Compton pho- 

tons constitute a distinctive bumplike feature in the dif- 
fuse background at ^300 MeV. For a boost factor near 
the middle of the allowed range (~1000), this bump will 
be only a factor of a few below the flux measured by 
EGRET. If the PAMELA excess is instead the result of 
the 300 GeV dark matter particle annihilating to muons, 
the resulting diffuse background will be less than ~ 10% 
of that measured by EGRET. 
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FIG. 2: The same as shown in Fig. [T] but after increasing the annihilation cross section (or boost factor) by the approximate 
maximum value consistent with the measurements of COMPTEL [ij and EGRET [l^l. For dark matter annihilating to 
electrons or tans, the inverse Compton photons provide a far more stringent constraint than prompt emission. 



If in addition to the PAMELA positron fraction, one 
would like to generate the electron spectrum reported by 
FGST with dark matter annihilations, we must consider 
a dark matter particle with a ~l-2 TeV mass and which 
annihilates largely to To normalize the annihila- 

tion rate to these signals, we must also introduce boost 
factors ranging from approximately 400-500 for a 1 TeV 
mass, to 1000 or more for a mass of 2 TeV [l^. Com- 
paring this to our results, we see that this would lead 
to a diffuse background constituting about one-fourth of 
the flux measured by EGRET around ~ 1 GeV. Simi- 
lar conclusions are reached if we consider dark matter 
that annihilates to muons through a light intermediate 



state 

With its superior angular resolution and greater expo- 
sure, FGST will study the diffuse extragalactic gamma- 
ray background between 100 MeV and hundreds of GeV 
with far greater precision than EGRET. In particular, 
FGST is expected to resolve and identify many of the 
individual sources that contribute to the diffuse isotropic 
background measured by EGRET. It has often suggested 
that a large fraction of the extragalactic diffuse back- 
ground measured by EGRET could consist of emission 
from blazars. As FGST is anticipated to resolve on the 
order of 10^ individual blazars [26] (EGRET, in contrast, 
only accumulated a catalog of 66 high-confidence blazars 
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FIG. 3: The maximum value of the boost factor (or equiv- 
alently, the maximum enhancement of the annihilation cross 
section relative to {(tv) — 3x 10~^^ cm'^/s) consistent with the 
diffuse isotropic background measured by the COMPTEL [ll| 
and EGRET 12] experiments. These constraints are espe- 
cially stringent for dark matter that annihilates directly to 
electrons or muons (for WIMP masses above 500 GeV in the 
latter case). 



over the duration of its mission [27|), the isotropic diffuse 
background that will be measured by FGST could be sig- 
nificantly smaller than that measured by EGRET j28| . 
If a significant fraction of the diffuse flux observed by 
EGRET originates from dark matter annihilations, then 
FGST may be capable of identifying this signal, espe- 
cially if a distinctive bumplike inverse Compton feature 
is present. Alternatively, if FGST resolves a large frac- 
tion of the background measured by EGRET, it will be 
able to place stringent constraints on the dark matter 
annihilation cross section to electrons or muons. 



SUMMARY AND CONCLUSIONS 

In this article, we have calculated the contribution to 
the extragalactic, isotropic gamma-ray background from 
annihilating dark matter particles, including the inverse 
Compton photons resulting from energetic electrons scat- 
tering with the cosmic microwave background. 

Although dark matter particles with typical thermal 
annihilation cross sections {{av) ^ 3 x 10~^^ cm^/s) are 
predicted to produce only a relatively small fraction of 
the isotropic gamma-ray and x-ray backgrounds observed 
by EGRET and COMPTEL, dark matter particles with 
larger cross sections (or equivalently, boost factors) could 
potentially be responsible for much or most of this flux. 



Dark matter scenarios capable of generating the cosmic 
ray signals observed by PAMELA, ATIC, or the Fermi 
gamma ray space telescope (FGST) predict dark matter 
particles which annihilate at a very high rate, and pri- 
marily to charged leptons. In such models, the gamma- 
ray background from inverse Compton scattering can be 
very significant, and is constrained by the measurements 
of COMPTEL and EGRET. 

If FGST resolves a significant fraction of the back- 
ground measured by EGRET into individual sources, it 
will be capable of placing very stringent limits on dark 
matter candidates with a high annihilation rate to lep- 
tons. If the dark matter annihilates directly to e~^e~, 
FGST may also be able to identify the distinctive inverse 
Compton peak at - 300 MeV x (TOx/300GeV)2 in 
the diffuse extragalactic gamma-ray spectrum. 

While we were in the final stages of completing this 
paper, Ref. appeared on the LANL preprint archive. 
Considering that they used the Moore density profile for 
dark matter halos, which accounts for approximately an 
order of magnitude flux enhancement over NEW profile, 
our results are in qualitative agreement. Once the halo 
profile is taken into account, Ref. |29| finds a flux that is 
a factor of a few larger than what we have found in our 
calculation. The possible sources of this discrepancy may 
include the multiplicity function of halos, the choice of erg 
and Us and the parametrization of the optical depth [g^ . 
We also note that our results for the prompt component 
of the signal are in very good quantitative agreement with 
those of Ref. 0] . This work has been supported by the US 
Department of Energy, including Grant No. DE-FG02- 
95ER40896, and by NASA Grant No. NAG5-f0842. 



APPENDIX 

Relativistic electrons transfer energy to photons in a 
process known as inverse Compton scattering. The rate 
of scattering of a relativistic electron of energy = 
jrUeC^ in a gas of monoencrgctic photons of energy i?rad 
and density A'lad to photons of energy Eic per energy 
interval dEic is 
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where e — ^^^^ ^ and (Tt is the Thompson cross section 
psj . In particular, electrons can scatter off of photons 
in the cosmic microwave background (CMB), which in 
our case might be considered nearly monoenergetic with 
energy density Uqmb = A^cmb£^cmb ~ 0.25eV/cm^ at 
present. The average energy of the scattered photon is 



F - ^(^^ \ F 
-C/IC — o ^CMB 

3 V me 



6 



3.4 X 10^2 GcV (1 



100 GcV 



(7) 



In order to obtain the spectrum of photons produced by 
a cascade of scatterings of an electron with initial energy 
Ec, we need to average over realizations. Using AE^ — 
Eic — EcMB oc E^ one can show that oc E^^^^^. 

Assuming the properties of the photon gas do not change 
while any given electron is losing most of its energy, the 
rate of energy transfer is 
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As a result of this process, a 10 GeV electron at redshift 
z = 1 will lose 90% of its energy to the CMB within about 
80 million years. Thus the vast majority of the energy 
that is deposited via dark matter annihilations into elec- 
trons gets almost immediately transferred into the cos- 
mic microwave background (relative to other timescales 
in the problem). At higher redshifts, or for higher energy 
electrons, the transfer of energy is even more rapid. 

Generally speaking while an electron/positron is losing 
its energy the CMB photons are redshifted. The charac- 
teristic time of energy loss due to IC scattering for a 100 
GeV electron/positron is tjc = 



TH 



miTTt « 3.9 X 10-s « 
4.35 X lO^^s, which is much less than the Hub- 



ble time, allowing us to neglect the redshifting of CMB 
photons. The spectral density of photons produced by 
an electron/positron of a given energy -^^{Eic,Ee) is 
then calculated by the Monte Carlo method. 

To obtain the spectral density of inverse Compton pho- 
tons for a given redshift per annihilation, one convolves 
the electron/positron density spectral density per anni- 
hilation with the spectral density of IC photons: 



dN 
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dE,^{Eic.E,)^m (9) 



The photon spectral density is then normalized to 
carry the same energy density as the electrons/positrons 
per annihilation. 
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